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Abstract 


While  Biot  theory  can  successfully  account  for  the  dispersion  observed 
in  sand  sediments,  the  attenuation  at  high  frequencies  has  been  ob¬ 
served  to  increase  more  rapidly  than  Biot  theory  would  predict.  In  an 
effort  to  account  for  this  additional  loss,  perturbation  theory  is  applied 
to  Biot’s  poroelastic  equations  to  model  the  loss  due  to  the  scattering 
of  energy  from  heterogeneities  in  the  sediment.  A  general  theory  for 
propagation  loss  is  developed  and  applied  to  a  medium  with  a  ran¬ 
domly  varying  frame  bulk  modulus.  The  theory  predicts  that  these 
heterogeneities  produce  an  overall  softening  of  the  medium  as  well  as 
scattering  of  energy  from  the  mean  fast  compressional  wave  into  inco¬ 
herent  fast  and  slow  compressional  waves.  This  theory  is  applied  to  two 
poroelastic  media:  a  weakly  consolidated  sand  sediment  and  a  consol¬ 
idated  sintered  glass  bead  pack.  The  random  variations  in  the  frame 
modulus  do  not  have  signihcant  effects  on  the  propagation  through 
the  sand  sediment  but  do  play  an  important  role  in  the  propagation 
through  the  consolidated  medium. 

BAGS  numbers:  43.30.Ma,  43.20.Gp,  43. 20. Jr 
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I.  INTRODUCTION 


For  sand  sediments  there  has  been  a  long-running  debate  about  the  physical  mechanisms 
by  which  sound  is  attenuated  as  it  passes  through  the  sediment.  This  debate  has  typically 
centered  on  two  explanations:  that  the  loss  is  due  to  viscous  flow  through  the  pores  formed 
by  the  sand  grains  which  compose  the  sediment  or  that  the  loss  is  due  to  the  interactions 
between  the  sand  grains  at  the  grain  contacts.  The  first  of  these  mechanisms  was  applied 
by  Stolh  using  the  theory  for  fluid-saturated  poroelastic  solids  developed  by  Biot^’^.  This 
theory  predicts  that  for  low  frequencies,  the  attenuation  will  increase  as  the  square  of  the 
frequency  while  at  higher  frequencies  it  will  increase  as  the  square  root  of  the  frequency. 
The  theory  also  suggests  that  the  change  in  attenuation  is  accompanied  by  an  increase 
in  sound  speed.  The  grain  contact  theory  was  proposed  by  Hamilton"^  and  later  revisited 
by  Buckingham^.  This  theory  predicts  that  the  attenuation  is  linear  in  frequency  over  all 
frequencies  while  the  sound  speed  exhibits  little  or  no  dispersion. 

A  number  of  experiments  have  been  performed  in  both  laboratory  and  ocean  sediments 
to  test  the  predictions  of  these  theories.  Hamilton  collected  a  number  of  the  measurements 
and  found  that,  in  the  aggregate,  the  attenuation  is  linear  in  frequency^.  This  observation 
has  been  supported  by  the  relatively  recent  sediment  acoustic  experiments  in  1999  and  2004 
(SAX99  and  SAX04  respectively)  for  high  frequencies®  as  well  as  recent  laboratory  experi¬ 
ments  examining  propagation  through  glass  bead  sediments^.  In  the  laboratory  experiments, 
while  Biot  theory  was  unable  to  capture  the  frequency  dependence  of  the  attenuation,  it  did 
predict  the  overall  scaling  of  the  attenuation  as  the  viscosity  of  the  pore  fluid  was  increased. 
Also  the  sound  speeds  measured  during  SAX99  and  SAX04  both  exhibit  dispersion  that 
for  the  mid-  to  high-frequencies  is  consistent  with  the  predictions  of  Biot  theory.  For  low 
frequencies,  there  remains  some  uncertainty  as  to  both  the  asymptotic  value  for  the  sound 
speed  as  the  frequency  goes  to  zero  and  the  frequency  dependence  of  the  attenuation®. 
Measurements  in  ocean  sediments  have  shown  that  for  low  frequencies  the  attenuation  is 
proportional  to  f"',  where  /  is  the  frequency  and  n  varies  from  1.5  to  2.0®. 


3 


These  measurements  have  led  to  attempts  to  create  modihcations  to  each  theory  to 
capture  aspects  of  the  other.  Buckingham’s  grain-to-grain  shearing  model  was  recently 
modified  to  incorporate  a  linear  loss  mechanism  at  the  grain  contacts  in  addition  to  the 
original  strain-hardening  mechanism®.  This  produces  a  linear  frequency  dependance  for  the 
attenuation  at  high-frequencies  and  a  dependence  similar  to  Biot  theory  in  the  mid-  to 
low-frequencies.  Likewise,  several  researchers  have  proposed  a  linear  loss  mechanism  at  the 
grain  contacts  which  would  manifest  as  complex,  frequency-dependent  frame  moduli  in  Biot 
theory  thus  creating  a  hybrid  of  the  two  approaches'^’^®. 

In  this  paper,  we  develop  an  alternative  theory  which  treats  the  sediment  as  a  poroelastic 
medium  using  Biot  theory  but  incorporates  losses  due  to  scattering  from  heterogeneities  in 
the  structure  of  the  sand  sediment.  These  heterogeneities  can  exist  in  any  of  the  parameters 
on  which  Biot  theory  depends  such  as  the  porosity  or  permeability.  For  a  fast  compressional 
wave  propagating  through  the  sediment,  energy  will  be  scattered  from  the  coherent  wave 
into  incoherent  fast  compressional  waves  when  the  wavelength  of  the  wave  is  on  order  of 
the  size  of  the  heterogeneity.  For  variations  on  the  order  of  a  few  grains,  these  losses  will 
occur  at  much  higher  frequencies  than  has  been  observed  in  the  experiments  discussed  above. 
However,  both  the  slow  compressional  wave  and  the  shear  wave  are  considerably  slower  than 
the  fast  wave  and  consequently  have  shorter  wavelengths.  When  these  wavelengths  are  on 
the  order  of  the  heterogeneities,  mode  conversion  can  occur  scattering  energy  from  the  fast 
compressional  wave  into  incoherent  slow  or  shear  waves  leading  to  additional  attenuation. 

This  type  of  attenuation  mechanism  has  been  explored  in  the  context  of  seismic  prop¬ 
agation  by  a  number  of  researchers.  Pride  and  Berryman  have  developed  a  model  for  the 
special  case  of  a  “double  porosity”  material  such  as  a  mixture  of  sand  and  clay  where  the 
heterogeneity  has  a  single,  dominant  length  scale^^’^^.  Pride  and  Masson  conducted  a  series 
of  numerical  simulations  to  measure  the  compressibility  of  a  heterogenous  poroelastic  mate¬ 
rial  which  has  a  self-affine  fractal  structure^®’^"^.  They  found  that  the  frequency  dependence 
of  the  attenuation  is  a  function  of  the  Hurst  exponent  which  characterizes  the  fractal  het¬ 
erogeneity.  Muller  and  Gurevich^^’^®  used  statistical  smoothing  of  the  Biot  Green’s  tensor  to 
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model  propagation  through  a  heterogenous  Biot  medium  for  frequencies  below  the  critical 


Biot  frequency, 


/3v 

u  <^UJb  = 

i^Pf 


(1) 


where  u  =  271  f  ^  (3  is  the  porosity,  rj  is  the  viscosity  of  the  pore  fluid,  p/  is  the  density  of 
the  pore  fluid,  and  n  is  permeability.  In  this  paper,  we  use  perturbation  theory  to  develop  a 
poroelastic  wave  equation  which  contains  the  hrst-order  scattering  effects  that  can  be  applied 
to  heterogenous  media  above  and  below  the  critical  frequency.  The  perturbation  technique 
used  was  hrst  developed  by  Karal  and  Keller  to  examine  scalar,  elastic,  and  electromagnetic 
propagation  through  random  media^^. 

Evidence  for  the  existence  of  heterogeneities  in  sand  sediments  can  be  found  in  recent 
work  in  the  physics  of  granular  materials.  Experiments^®  and  simulations^®  have  shown  that 
the  complexity  of  a  random  packing  of  sand  grains  or  beads  leads  to  the  formation  of  force 
networks  or  “chains”  where  a  small  subset  of  the  grains  will  carry  a  majority  of  any  load 
applied  to  the  medium  as  well  as  the  overburden  pressure  of  the  grains  themselves.  Hertz- 
Mindlin  theory  relates  the  forces  at  the  grain  contacts  to  the  frame  moduli  of  the  medium^® 
and  one  would  expect  that  in  applying  that  theory,  these  force  chains  should  manifest  as 
small  scale  heterogeneities  in  frame  moduli.  Likewise,  recent  analysis  of  simulations  of 
random  packings  of  spheres^^  as  well  as  measurements  on  sand  sediments^^  indicate  that 
the  porosity  should  also  exhibit  random  spatial  variations.  In  this  paper,  we  will  only 
consider  variations  in  the  frame  bulk  modulus  to  simplify  the  development  of  the  theoretical 
framework  which  will  be  used  to  consider  variations  in  the  porosity  in  a  subsequent  paper. 

Note  that  while  there  has  been  some  debate  over  whether  or  not  the  Biot  slow  wave 
can  be  observed  in  a  sand  sediment^®^^®,  the  theory  presented  here  is  not  predicated  on 
the  existence  of  a  propagating  Biot  slow  wave.  For  small  values  of  the  frame  moduli,  Biot 
theory  predicts  that  the  Biot  slow  wave  should  be  highly  attenuated  making  it  difficult  if 
not  impossible  to  observe^®.  However,  it  should  still  be  possible  to  scatter  energy  into  this 
type  of  wave,  in  which  case,  that  energy  will  be  quickly  converted  into  heat. 

In  Section  II,  we  begin  by  expressing  the  heterogenous  poroelastic  equations  in  a  form 
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to  which  we  can  apply  perturbation  theory,  the  details  of  which  are  developed  in  Section  III. 
The  general  properties  of  the  resulting  perturbed  wave  equation  including  the  limits  of 
applicability  and  causality  are  examined  in  Section  IV.  In  Section  V,  the  perturbation  theory 
is  applied  to  the  SAX99  and  sintered  glass  bead  sediments  and  the  results  are  discussed  in 
Section  VI. 


II.  HETEROGENOUS  BIOT’S  POROELASTIC  EQUATIONS 

If  the  frame  moduli  and  the  porosity  are  allowed  to  vary  as  functions  of  position  within 
a  poroelastic  medium,  Biot’s  poroelastic  equations  can  be  expressed  as^^ 


ai[(R-2/ife)V-u-C'V-w]  + 

dj  [nb  {diUj  +  djUi)]  =  -u'^pui  +  u^pfWi 


(2) 


and 


di  [UV  ■  u  —  MV  •  w]  =  —u'^pfUi  +  uj^p*Wi, 


(3) 


where  pb  is  the  shear  modulus  of  the  sediment  frame,  u  is  the  displacement  field  of  the 
frame,  w  =  /?  (u  —  Uf)  is  the  relative  displacement  of  the  pore  fluid  to  the  frame,  Uf  is  the 
displacement  field  of  the  pore  fluid,  Pg  is  the  sediment  particle  mass  density,  and 


D-Kb  3  ’ 


C  = 


Kg  {Kg  -  Kb) 

D-Kb 


K^ 

M=  " 


D  =  Kr. 


D-Kb' 

-’(S- 


and 


P  =  -  P)  Pg  +  Ppf, 


KLO 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 
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where  Kg  is  the  bulk  modulus  of  the  individual  grains,  is  the  bulk  modulus  of  the  sediment 
frame,  Kf  is  the  bulk  modulus  of  the  pore  fluid,  a  is  the  tortuosity,  and  the  parameter  F 
represents  the  deviation  from  Poiseuille  flow  as  frequency  increases.  The  expression  for  F  is 
given  and  discussed  in  Ref.28. 

The  two  equations  can  be  consolidated  into  a  single,  matrix  equation  as 


V|([K] 


where 


2[p])V-U]  +  V  X  ([p]  V  X  U)  + 
2  [il]  V^U  +  2  (V  [p]  ■  V)  U  + 

[p]  = 


[K] 

[P] 


H 

C 

P 

Pf 


0 

-c 

-M 

-Pf 

-P* 


fXb  0 
0  0 


(10) 

(11) 

(12) 

(13) 


and 


(14) 


which  for  the  homogenous  case,  when  [K]  =  [Kq]  and  [p]  =  [po]  are  constants,  rednces  to 


([Ko]  -  [po])  V  (V  ■  U)  +  [po]  V^U  +  [po]  =  0.  (15) 


The  homogenous  equation  can  be  solved  in  terms  of  potentials  to  obtain  the  fast  and  slow 
compressional  wavennmbers  and  the  shear  wavenumber. 

In  the  following,  we  are  going  to  look  at  heterogeneities  due  to  the  frame  bulk  mod¬ 
ulus  only.  For  this  case,  the  frame  shear  modulus  and  the  porosity  are  constant  and  the 
poroelastic  eqnation  becomes 

V  ([K]  V  ■  U)  -  [po]  V  X  (V  X  U)  +  [po]  cu'U  =  0.  (16) 
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III.  PERTURBATION  THEORY 


In  this  paper  we  consider  a  poroelastic  medium  in  which  the  frame  bulk  modulus  varies 
randomly  with  position, 


/V(r)  =  A'j„  +  Vl',(r), 


(17) 


where  Kbo  =  {Kb  (r))  is  the  average  over  an  ensemble  of  realizations  of  the  medium  (denoted 
by  (■)),  SKb(r)  is  the  local  fluctuation  of  the  bulk  modulus,  and  {SKb{r))  =  0.  The  bulk 
modulus  enters  the  heterogenous  poroelastic  equations  through  the  Biot  coefficients  in  the 
matrix  [K],  To  account  for  the  random  variations,  we  will  expand  this  matrix  as  a  Taylor 
series  about  the  mean  value. 


1K|  =  [Ko]  +  |Fbi]  SIQ  +  -  |Fb2]  Vlt  +  ■  ■  ■ . 


(18) 


where 


Fbn 


<9"[Kc 

dKJ! 


(19) 


This  expansion  can  now  be  substituted  into  the  poroelastic  equation  given  by  Eq.  (16).  If 
the  heterogenous  poroelastic  equation  is  rewritten  as 


£([K]).U  =  0, 


(20) 


substitution  of  the  expanded  expressions  yields, 

(To  +  Ti  +  T2  +  •••)■  U  =  0,  (21) 

where  To  ■  U  =  0  is  the  homogenous  poroelastic  equation  given  by  Eq.  (15).  We  will  also 
expand  the  displacement,  U,  as 

U  =  (U)  +  U„  (22) 

where  (U)  is  the  mean  coherent  held  and  Ug  is  the  incoherent,  scattered  held  for  which 
(Ug)  =  0.  Substitution  of  Eq.  (22)  into  Eq.  (21)  yields. 

To  ■  (U)  +  To  ■  Ug  +  Ti  ■  (U)  +  Ti  ■  Ug 


+  T2  ■  (U)  =  0, 


(23) 


where  we  have  neglected  all  terms  higher  than  second  order. 

Our  goal  is  to  obtain  a  solution  for  the  coherent  held,  (U).  Since  (£1)  =  0,  the  mean  of 
Eq.  (23)  yields 

£0  ■  (u)  +  (£1  ■  Us)  +  (£2)  ■  (U)  =  0.  (24) 

Subtracting  Eq.  (24)  from  Eq.  (23)  and  neglecting  all  terms  higher  than  hrst  order,  we  hnd 
the  equation  for  Ug  in  terms  of  (U), 

£0  ■  Us  =  -£i  ■  (U) .  (25) 

This  equation  can  be  solved  using  the  poroelastic  Green’s  tensor,  which  is  a  solution  to  (see 
Appendix  A) 

£(r)  ■  0(r,r')  = -5(r,r')  [I]  J,  (26) 

to  yield  the  solution  for  Us, 

Us  {r)  =  (r,  r')  ■  £1  (r')  ■  (U  (r'))  dh'.  (27) 

Substituting  this  expression  back  into  Eq.  (24),  we  hnd  the  hnal  equation  for  a  poroelastic 
wave  propagating  through  a  heterogenous  medium, 

£o(r)  ■  {U)(r) 

+  Ai  (r)  ■  0  (r,  r')  ■  £,  (r')  ■  (U  (r'))  <i’rA 

+  (£.(r)>.{U(r))  =  0.  (28) 

This  result  is  identical  to  Karal  and  Keller’s  Eq.  (14)  in  Ref.  17.  for  the  general  case  of 
scalar,  elastic,  or  electromagnetic  wave  propagation  in  a  random  medium. 

In  order  to  apply  Eq.  (28)  to  the  case  of  poroelastic  waves  propagating  through  a 
heterogenous  medium,  we  begin  by  substituting  Eq.  (18)  into  the  heterogenous  poroelastic 
equation,  Eq.  (16).  The  result  can  be  separated  into  three  parts:  £o  ■  U,  which  is  identical 
to  Eq.  (16)  with  the  replacement  of  [K]  by  [Kq], 

£i  ■  U  =  [Fbi]  V  {6KbV  ■  U)  (29) 
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and 


=  \  [Fb2]  V  {6K^,V  ■  U)  .  (30) 

The  last  term  in  Eq.  (28)  can  be  determined  by  substituting  U  =  (U)  in  Eq.  (30)  and 
taking  the  mean  of  the  result.  Since  (SKh)  =  0,  this  yields, 

{£2(r)>.{U)  =  ilFt2|ff^V(V.(U)),  (31) 

where  =  (SK^)  is  the  variance  of  the  perturbation  in  the  frame  bulk  modulus. 

Consider  the  integral  factor  in  Eq.  (28)  involving  the  Green’s  tensor.  This  factor  can  be 
viewed  as  the  field  due  to  a  source  £i  ■  U.  Eq.  (29)  shows  that  this  source  is  the  gradient 
of  a  scalar  held,  from  which  Eq.  (A15)  leads  to 

j  0(r,rO-i:i(r')-(U(r'))dV  (32) 

^  r 

=  [Fbi]  klV  /  gn  (R)  SK,  (rO  V'  ■  (U  (r'))  dh'. 

n=l  '' 

where  V'  denotes  differentiation  with  respect  to  r'  and  the  Qn  {R)  are  dehned  in  Appendix  A. 
Note  that  as  a  consequence  of  the  source  being  the  gradient  of  a  scalar  held,  there  is  no  mode 
scattering  into  shear  waves.  Operating  on  this  expression  with  Hi  and  taking  the  mean  will 
yield  the  second  term  in  Eq.  (28), 

(r)  ■  y  ^  (r,  r)  ■  Hi  (E)  ■  (U  (r'))  = 

2 

-  IFbil  Y.  1= J  |F«I 

n=l 

|y,V  Ulh  (r)  j  9„  (R)  SIQ  (r')  V'  ■  (U  (r')) 

+  <T,\V(V.{U(r)))}  (33) 

where  we  have  used, 

V^gn{R)  =  -klgn{R)-S{R).  (34) 

In  Eq.  (33),  the  mean  is  taken  only  over  the  product  of  the  frame  modulus  varia¬ 
tions  and  the  resulting  expression  is  the  covariance  of  the  modulus  variations,  Cbb  (r,  rO  = 
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{SKh  (r)  SKf,  (r')).  If  we  assume  that  the  variations  in  the  modulus  are  statistically  homoge¬ 
nous  (spatially  stationary)  and  isotropic,  then  the  covariance  should  depend  only  on  the 
distance  R  =  |r  —  r'|  and  Cbb  (r,  r')  =  Cbb  (R)- 

In  order  to  solve  the  perturbation  equation,  we  will  assume  a  plane  solution  given  by 

(U(r))  =  (35) 

With  this  expression  for  the  mean  held  and  using  Eqs.  (15),  (31),  and  (33),  Eq.  (16) 
becomes 

-([Ko]  -  [p-o])  k  (k  •  A)  -  [po]  A  [po]  cn^A 
2 

+  [Fbil5;[G„l  MfcJk(k.A)  X 

n=l 

I  (y,9„  (R)  +  HR))  Cu.  (R)  e-" 

-l|Fb2l<T?tk(k.A)=0,  (36) 

where  we  have  substituted  R  =  r  —  r'  and  integrated  with  respect  to  R.  Evaluating  the 
integral  of  the  delta  function  and  using  the  relation  given  in  Eq.  (A16),  we  can  rearrange 
the  terms  to  express  the  wave  equation  in  the  much  simpler  form, 

-  ([Keff]  -  [p-o])  k  (k  ■  A)  -  [po]  A  -P  [po]  cn^A  =  0,  (37) 

where  the  effects  of  the  bulk  frame  modulus  variations  on  the  wave  propagation  have  been 
folded  into  the  effective  Biot  coefficient  matrix, 

[Keff]  =  [Ko] 

2 

-  IFbi]  ^  |G„]  |F„]  k‘I„ 

n=l 

-  (|Fbil[K„l-MFbil-l[Fb2l])<  (38) 

where 

In  =  j  9n  {R)  Cbb  {R)  (39) 


11 


When  the  variance  of  the  modulus  variations  goes  to  zero,  the  effective  matrix  reduces  to  [Kq] 
as  one  would  expect.  The  second  term  which  contains  the  integration  over  the  covariance 
will  lead  to  dispersion  and  an  increase  in  attenuation  as  a  function  of  frequency  as  will 
be  seen  in  the  examples  shown  in  Section  V.  The  last  term  in  this  expression  does  not 
depend  on  the  form  of  the  covariance  nor  does  it  have  any  frequency  dependence.  This  term 
represents  the  softening  of  the  medium  due  to  the  variations  in  the  frame  modulus  which  is 
discussed  in  more  detail  in  the  next  section  and  in  Appendix  B. 

The  result  in  Eq.  (37)  has  three  solutions:  a  shear  wave,  a  slow  compressional  wave,  and 
a  fast  compressional  wave.  For  the  shear  wave,  the  displacement  of  the  wave  is  perpendicular 
to  the  wave  vector,  k,  and  k  (k  •  A)  =0.  The  random  variations  in  the  bulk  frame  modulus 
have  no  effect  on  the  shear  wave  and  the  resulting  expression, 

-/c^  [po]  A  +  [po]  cu^A  =  0,  (40) 


is  identical  to  the  result  for  the  shear  wave  in  a  homogenous  medium.  For  the  compressional 
waves,  the  displacement  is  parallel  to  the  wave  vector  and  Eq.  (37)  reduces  to 

—  [Kgff]  /c^A  +  [po]  ci;^A  =  0,  (41) 


which  has  two  solutions  since  it  is  a  2  x  2  matrix  equation  (Note  that  the  equation  for  the 
shear  wave  has  only  a  single  solution  since  only  the  hrst  element  in  [po]  is  non-zero).  For 
the  remainder  of  this  paper,  we  will  focus  on  the  fast  compressional  wave  solution,  ki. 

The  integral  in  [Keff]  given  by  Eq.  (39)  can  be  partially  solved  by  separating  it  into  its 
radial  and  surface  components. 


In  —  j  Qn  {R)  Cbb  (i?)  /  e 


— ik-R 


dSdR. 


(42) 


The  integration  with  respect  to  dS  can  be  performed  yielding, 

f  e-‘'‘-^dS  ^^sinkR. 


(43) 


Substituting  this  result  back  into  Eq.  (42)  and  substituting  the  expression  in  Eq.  (A7)  for 
the  scalar  green’s  function,  gn  (R),  the  hnal  integral  becomes. 


k 


roo 


jikn  R 


Cbb  (R)  sin  (kR)  dR. 


(44) 
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The  fast  compressional  wave  propagating  through  a  poroelastic  medium  with  a  randomly 
varying  bulk  frame  modulus  can  now  be  determined  from  Eq.  (41)  where  [Keff]  is  given  by 
Eq.  (38)  and  is  given  by  Eq.  (44). 

In  the  low  frequency  limit,  when  lo  satishes  Eq.  (1),  the  effective  wavenumber  determined 
from  Eq.  (41)  has  a  much  simpler  form.  In  this  limit,  p*  — >  irf/nuj^  and  the  unperturbed 
compressional  wavenumbers  become 


ki  =  u 


and 


(45) 


where  N  =  M  —  / H  and  ki  is  real.  The  slow  wave,  /c2,  in  this  low  frequency  limit,  is 
a  solution  to  a  diffusion  equation  (See  Eq.  (5a)  and  (5b)  in  Ref.  26.)  In  this  limit,  if 
we  consider  only  scattering  from  fast  into  slow  waves  and  assume  that  the  compressional 
wavenumbers  satisfy 

<  1, 

then  the  effective  fast  compressional  wavenumber  reduces  to 


where 


and 


keff  ui 


H, 


eff 


Heff  =  H-  A2ai  -  A,kl  I  Ctb  (R) 


_  M 

dKh  oKh 


A2  —  Ai  — 


dH 

m,, 

Id'^H 


H 


dC 

Ml 


2  dK^ ' 


(47) 

(48) 

(49) 

(50) 


This  result  is  equivalent  to  the  effective  wavenumber  in  Eqs.  (49),  (50),  and  (51)  of  Ref.  16 
by  Miiller  and  Gurevich  when  there  are  no  variations  in  the  shear  modulus. 


IV.  GENERAL  PROPERTIES  OF  COMPRESSIONAL  SOUND  SPEED 
AND  ATTENUATION 

We  will  explore  the  behavior  of  the  perturbed  fast  compressional  wave  for  two  types 
of  media:  the  SAX99  sand  sediment  and  a  water-saturated  sintered  glass  bead  pack.  The 
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1.00 


(b) 

FIG.  1.  (a)  Fast  compressional  sound  speed  and  (b)  attenuation  determined  from  Eq.  (41) 
with  [Kgff]  given  by  Eq.  (54)  as  a  function  of  the  variance,  of  the  moduli  variations 
normalized  by  K^.  Both  the  sound  speed  and  the  attenuation  have  been  normalized  by 
their  respective  values  when  Chh  =  0.  These  results  were  calculated  for  /  =  1  kHz  using 
the  SAX99  sediment  parameters  in  Table  I  while  the  value  of  the  mean  bulk  modulus  was 
increased  from  4.36  x  10^  Pa  for  the  solid  line,  to  4.36  x  10®  Pa  for  the  short  dashed  line, 
and  to  2.18  x  10®  Pa  for  the  long  dashed  line. 
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SAX99  sediment  is  a  well-characterized  sand  sediment  which,  like  most  sand  sediments,  has 
a  small  frame  bulk  modulus.  The  values  for  the  SAX99  sediment  given  in  Table  I  were  taken 
from  Table  I  of  Ref.  6.  The  sintered  bead  pack  was  chosen  since  it  has  a  much  larger  frame 
bulk  modulus  and  has  been  studied  extensively  as  a  poroelastic  medium^®’^®’^°.  A  majority 
of  the  material  parameters  for  the  sintered  glass  bead  pack  were  taken  from  Sample  2  in 
Table  III  of  Ref.  30  while  the  fluid  parameters  were  chosen  to  correspond  to  the  water 
parameters  measured  during  SAX99.  These  two  media  allow  us  to  understand  the  effects  of 
scattering  from  heterogeneities  for  both  a  weak  and  a  strong  frame. 


A.  Softening  of  the  poroelastic  medinm 

In  order  to  understand  the  behavior  of  the  fast  compressional  wave  in  a  heterogenous 
medium,  we  begin  by  considering  the  solution  to  the  perturbed  wave  equation  for  a  poroe¬ 
lastic  medium  where  the  bulk  frame  modulus  variations  are  related  by  the  exponential 
correlation  function, 

=  (51) 


where  L  is  the  characteristic  length  scale  of  the  bulk  modulus  variations.  The  three- 
dimensional  Fourier  transform  of  this  correlation  function  yields  the  commonly  encountered 
power  spectrum. 


Wbb  {k) 


7r2L(A;2  +  L-2)2‘ 


(52) 


The  integral  given  in  Eq.  (44)  can  be  solved  analytically  for  the  exponential  correlation 
function  to  yield. 


In 


(j 


2 

bb 


ikn)^  +  k'^ 


(53) 


In  the  following  section,  we  will  explore  this  result  in  detail,  but  here  we  will  consider  the 
limiting  case  as  L  — 0.  This  corresponds  to  a  poroelastic  medium  in  which  the  frame 
modulus  is  completely  random  between  any  two  points  in  the  medium.  In  this  case,  /n  =  0 
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FIG.  2.  Values  of  for  which  [Kgff]  remains  equal  to  [Kq]  for  the  SAX99  sediment 

parameters  (solid  line)  and  the  sintered  glass  bead  parameters  (dash-dot  line)  as  Kb  is 
increased  by  three  orders  of  magnitude. 

and  the  effective  Biot  coefficient  matrix  reduces  to 

[Keff]  =  [Kq]  -  ^[Fbl]  [Kq]  ^  [Fbl]  -  -  [Fb2]^  (54) 

As  shown  in  Appendix  B,  this  expression  for  [Keff]  is  the  effective  response  of  a  poroelastic 
medium  with  a  heterogenous  frame  modulus  under  static  compression.  The  elastic  response 
of  the  medium  is  dominated  by  the  weaker  elements  of  the  frame  and  consequently  the 
heterogenous  frame  is  softer  than  the  homogenous  medium. 

This  softening  has  two  important  consequences.  First  it  produces  a  decrease  in  the 
sound  speed  that  is  proportional  to  the  variance  as  seen  in  Figure  1(a).  The  sound  speed 
has  been  calculated  using  the  SAX99  sediment  parameters  given  in  Table  I  for  the  solid  line. 
The  bulk  frame  modulus  was  increased  by  an  order  of  magnitude  for  the  short-dashed  line 
and  two  orders  of  magnitude  for  the  long-dashed  line.  This  softening  of  the  medium  also 
leads  to  an  increase  in  the  attenuation  as  the  variance  increases.  The  attenuation  which 
corresponds  to  the  sound  speed  plotted  in  Figure  1(a)  is  shown  in  Figure  1(b). 

This  dependence  on  the  variance  for  both  the  sound  speed  and  attenuation  may  have  an 
important  influence  on  the  choice  of  frame  modulus  used  in  Biot  calculations.  For  a  static 
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measurement  of  the  elastic  response  of  a  sediment,  if  the  frame  modulus  is  heterogenous, 
the  value  determined  for  the  bulk  frame  modulus  will  not  be  the  mean  frame  modulus  but 
instead  a  smaller  value.  Likewise  for  a  measurement  of  the  sound  speed,  the  frame  modulus 
determined  from  a  best  £t  of  the  data  would  again  underestimate  the  mean  value  for  the 
frame  modulus.  If  K^q  is  the  mean  bulk  modulus  for  a  homogenous  sediment,  it  is  possible 
to  hnd  values  for  and  that  satisfy, 

[Kn]{K,,ai)  =  [Ko]{Kto,0),  (55) 

where  the  effective  Biot  coefficients  as  a  function  of  are  equal  to  the  coefficients  when 
=  0.  This  equation  can  be  solved  for 

2  {D  -  K,)  (K,  -  K,o)  f 

= - (D^) - <“> 

This  result  is  plotted  in  Figure  2  for  both  the  SAX99  and  sintered  glass  bead  sediments.  For 
the  SAX99  sediment,  as  the  bulk  modulus  is  increased,  the  normalized  variance  comes  very 
close  to  unity  and  then  decreases  while  for  the  sintered  glass  beads,  the  normalized  variances 
increases  to  almost  1.2  before  decreasing  to  zero.  The  variance  becomes  zero  when  K),  =  D, 
where  D  is  given  in  Eq.  (7).  At  this  point  H,C,  and  M  diverge  and  both  the  perturbed 
and  unperturbed  wavenumbers  equal  zero. 

For  both  media,  Eq.  (56)  predicts  that  it  is  possible  to  choose  values  of  and  Ki, 
where  is  almost  equal  to  or,  for  the  case  of  the  sintered  glass  beads,  larger  than  It 
is  worth  considering  at  this  point  whether  or  not  large  values  could  be  realistic  for  actual 
sediments  or  other  poroelastic  media.  The  frame  moduli  are  constrained  to  be  positive  by  the 
nonnegative  character  of  the  strain  energy  for  the  medium^.  If  the  probability  distribution 
for  the  modulus  variations  is  approximately  Gaussian,  a  standard  deviation,  abb,  on  the 
order  of  the  mean  value  for  the  frame  modulus  would  imply  that  at  some  locations  within 
the  medium,  the  local  frame  modulus  would  have  to  be  negative,  which  would  violate  the 
constraint  of  the  strain  energy.  However,  the  distribution  of  contact  forces  within  real 
and  simulated  granular  materials  have  been  found  to  be  approximately  exponential.  The 
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probability  distribution  for  the  frame  modulus  variations  may  exhibit  a  similar  kind  of 
distribution,  in  which  case,  the  values  for  in  Fig.  2  may  be  reasonable. 


B.  Limits  of  applicability 

These  large  values  for  the  variance  also  bring  into  question  the  validity  of  the  assumptions 
behind  the  derivation  in  Section  III,  specihcally  whether  it  is  sufficient  to  account  for  only 
the  second-order  scattering  contributions  to  the  mean  propagating  held.  As  the  variance 
increases,  it  may  be  necessary  to  incorporate  the  effects  of  fourth  order  contributions.  This 
can  be  done  formally  by  retaining  the  fourth  order  terms  in  Eq.  (23)  and  following  the  same 
procedure  used  in  Eq.  (24)- (28).  While  it  may  be  possible  to  solve  the  resulting  expression, 
the  contributions  of  these  higher  order  terms  can  be  estimated  in  a  more  straightforward 
way  to  yield  a  necessary  condition  which  must  be  met  for  the  second  order  theory  to  be 
valid. 

The  second  approach  to  testing  the  validity  of  the  second  order  approximation  follows 
the  approaches  of  Rytov^^  and  Miiller^®.  They  examine  the  next  approximation  for  the 
kernel  of  mass  operator,  which  replaces  the  unperturbed  Green’s  function  with  the  Green’s 
function  they  obtained  from  the  original  approximation.  For  the  current  theory,  it  can  be 
shown  that  an  inhnite  subset  of  terms  in  the  full  expansion  corresponding  to  Eq.  (23)  can 
be  solved  by  replacing  the  Green’s  tensor  in  Eq.  (38)  with  the  mean  Green’s  tensor,  0  (r,  r') 
for  the  heterogenous  medium.  The  mean  Green’s  tensor  can  be  found  by  considering  the 
inhomogenous  poroelastic  equation  with  a  spatially  varying  frame  modulus, 

V  ([K]  V  ■  0)  +  [po]  V  X  (V  X  0)  + 

[Po]o;20  =  -5(R)[|]  J.  (57) 

We  showed  in  Section  III  that  the  homogenous  wave  equation  for  the  mean  held  can  be 
expressed  in  terms  of  the  unperturbed  wave  equation  with  [Kq]  replaced  by  [Keff] .  It  follows 
that  the  inhomogenous  wave  equation  can  be  approximated  in  the  same  way  and  0  is  a 
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solution  to 


[Keff]  V  (V  ■  0)  +  [fio]  V  X  (V  X  0)  + 

[po]o;20  = -5  (/?)[!]  3.  (58) 

The  mean  Green’s  tensor  therefore  has  the  same  form  as  the  unperturbed  Green’s  tensor 
given  by  Eq.  (A3)  with  [Gi]  and  [G2]  replaced  by 

- ’ 

/il  n.2 

and 

H  - ’ 

where  ki  and  k2  are  the  mean  compressional  wavenumbers  determined  from  Eq.  (41).  The 
mean  compressional  wavenumbers  also  replace  the  unperturbed  wavenumbers  in  Eq.  (A7). 

The  next  approximation  for  the  fast  compressional  wave  can  be  found  by  substituting 
0  for  0  in  Eq.  (28).  The  solution  to  the  new  equation  proceeds  as  before  to  yield  the  new 
effective  Biot  coefficient  matrix, 

k2]  =  [K„| 

2 

-  [Uii  E  [SJ  IF-'I 

n=l 

-  (|Fbi]  |K.„|-‘  1F„]  -  i  4,  (61) 

which  can  be  substituted  into  the  compressional  wave  equation  to  get  the  new  effective 
wavenumber,  k[^k 

The  condition  which  must  be  met  is  that  the  contribution  of  this  higher  correction  to  the 
wavenumber  be  small  compared  to  the  second  order  correction  to  the  wavenumber,  which 
can  be  expressed  as 

-  ki 
ki  —  kij 

This  constraint  is  a  necessary  condition  which  must  be  met  in  order  for  the  second-order 
approximation  to  apply  to  heterogenous  poroelastic  medium. 
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C.  Causality 


For  the  low-frequency  limits  of  the  mean  compressional  wavenumber,  given  by  Eqs.  (47)- 
(50),  Muller  and  Gurevitch  showed  that  this  result  satished  a  twice-subtracted  dispersion 
relation  and  hence  was  causah®.  For  the  general  wavenumber  result,  which  is  applicable  at  all 
frequencies,  we  cannot  show  analytically  that  the  complex  wavenumber  satishes  a  dispersion 
relation.  It  is  possible  however  to  numerically  integrate  a  twice-subtracted  dispersion  relation 
for  the  real  and  imaginary  parts  of  the  mean  wavenumber  on  a  case  by  case  basis.  If  we 
dehne  an  index  of  refraction,  n{uj)  =  coki/uj,  where  cq  is  an  arbitrarily  chosen  real  sound 
speed,  then  the  twice  subtracted  dispersion  relation  that  must  be  evaluated  becomes 


R 


n  (u)  = 

2 

TT 


cut 


-|-  n  : 


R 


UJr, 


—  ool)  x 

“  {u')  uj'duj' 


+ 


(63) 


Jo  (cn'2  -  uf)  (cn'2  -  (a;'2  -  a;^)  ’ 

where  n  =  +  in^ .  The  indices  rii  =  n  (cui)  and  n2  =  n  {002)  are  evaluated  at  the  arbitrary 

frequencies  uii  and  UJ2.  Choosing  cq  =  Ui/ki  (cui),  this  equation  can  be  numerically  evaluated. 

For  the  examples  given  in  the  next  section,  the  dispersion  relation  given  above  was  found 
to  be  satished  to  within  the  numerical  errors  of  the  calculation.  This  does  not  prove  that 
the  perturbed  wavenumber  is  causal  in  all  cases,  but  rather  provides  a  strong  indication  that 
the  theory  is  causal. 


V.  SOUND  SPEED  AND  ATTENUATION  FOR  SPECIFIC  CORRELATION 
FUNCTIONS 

In  order  to  model  a  heterogeneous  medium  using  the  theory  developed  in  Section  III,  the 
correlation  function  for  the  bulk  modulus  variations  must  be  known  or  determined  from  a  ht 
to  existing  data.  In  this  section  we  consider  two  examples  of  correlation  functions  frequently 
used  in  acoustics  and  geophysics. 

The  hrst  of  the  correlation  functions  to  be  considered  is  the  exponential  correlation 
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(a) 


(b) 

FIG.  3.  Examples  of  (a)  sound  speed  ratio  and  (b)  attenuation  for  the  SAX99  sediment 
assuming  an  exponential  correlation  for  the  frame  bulk  modulus  variations.  The  correlation 
length  is  held  hxed  at  L  =  1  mm  while  the  variance  and  the  mean  bulk  modulus  are 
increased.  The  variance  is  shown  in  the  hgures  and  the  mean  bulk  modulus  is  iFf,  =  l.lSiFfeo 
(long-dashed  line),  iFf,  =  lAlK^o  (short-dashed  line),  and  Kb  =  2.04iFbo  (long-short  dashed 
line).  In  (a)  the  sound  speed  ratio  is  dehned  as  the  sound  speed  of  the  sediment  normalized 
by  the  sound  speed  of  the  pore  fluid.  In  (b)  the  lower  set  of  curves  is  the  contribution 
of  the  scattering  from  the  fast  compressional  waves  into  the  slow  compressional  waves  to 
the  attenuation  of  the  mean  fast  compressional  wave.  The  solid  line  in  each  figure  is  the 
prediction  of  Biot  theory  with  no  variation  in  the  bulk  modulus  and  Kb  =  Kbo- 
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function  which  was  introduced  in  the  previous  section.  Using  the  solution  of  the  integral  for 
the  exponential  correlation  given  in  Eq.  (53),  the  effective  Biot  coefficient  matrix  becomes, 


[Keff]  —  [Kq] 


C, 


bb 


IFbil  ^  |G„]  |F 


n=l 


bb 


:Fbi]  [K( 


1-1 


bl 


'bl 


K 


(L  1  -  ikn)  +  k'^ 


i  1F« 


(64) 


Note  that  the  wavenumber  that  we  seek  to  determine,  k,  appears  in  the  denominator  of  the 
second  term  in  [Kgff] .  In  order  to  solve  the  wave  equation,  we  approximate  the  wavenumber 
in  the  denominator  by  the  unperturbed  fast  wavenumber,  ki.  The  validity  of  this  approx¬ 
imation  can  be  evaluated  by  substituting  ki  into  the  denominator  and  iterating  until  the 
solution  converges.  For  the  examples  given  here,  this  approximation  yields  a  solution  that 
is  not  significantly  different  than  the  iterated  solution. 

Examples  of  the  sound  speed  and  attenuation  evaluated  using  Eq.  (64)  as  a  function  of 
are  shown  in  Figure  3  for  the  SAX99  sediment.  In  each  of  these  cases,  the  correlation  length 
was  held  fixed  at  L  =  1  mm  and  K),  was  increased  according  to  Eq.  (56)  to  insnre  that  the 
low  freqnency  sonnd  speed  matches  the  unperturbed  result.  For  the  values  of  and 
chosen,  the  sonnd  speed  and  the  attenuation  are  only  weakly  affected  by  the  scattering  from 
fast  into  slow  compressional  waves  while  the  fast  to  fast  wave  scattering  occurs  at  much 
higher  frequencies.  The  magnitnde  of  the  scattering  contribution  to  the  attenuation  can  be 
seen  in  Fig.  3(b)  where  the  set  of  cnrves  in  the  lower  right-hand  corner  of  the  figure  were 
determined  by  subtracting  the  nnpertnrbed  attennation  from  the  perturbed  result.  Energy 
begins  to  scatter  from  the  mean  fast  compressional  wave  into  the  slow  wave  when  k2L  ^  1 
and  then  the  attenuation  dne  to  this  mode  conversion  levels  off  at  k2L  k,  10. 

As  noted  in  Fig.  2,  it  is  possible  to  choose  larger  valnes  for  Kiy  and  that  satisfy  the 
relation  given  in  Eq.  (56).  However,  if  we  evalnate  the  condition  given  in  Eq.  (62)  which 
determines  whether  the  second  order  scattering  theory  is  snfficient  to  model  the  propagating 
field,  we  find  that  the  right  hand  side  of  Eq.  (62)  is  eqnal  to  0.1,  0.25,  and  0.65  for  the  three 
cases  examined  in  the  figure.  We  have,  therefore,  already  begun  to  violate  the  inequality  in 


22 


FIG.  4.  Examples  of  sound  speed  ratio  and  attenuation  for  the  sintered  glass  beads  as¬ 
suming  an  exponential  correlation  for  the  frame  bulk  modulus  variations.  In  (a)  and  (b), 
the  correlation  length  is  held  hxed  at  L  =  1.6  mm  while  the  variance  and  the  mean  bulk 
modulus  are  increased.  The  variance  is  shown  in  the  figures  and  the  mean  bulk  modulus  is 
Ki,  =  1.009i^bo  (long-dashed  line),  Ki,  =  1.021if'feo  (short-dashed  line),  and  Ki,  =  1.045i^bo 
(long-short  dashed  line).  In  (c)  and  (d),  the  variance  and  mean  bulk  modulus  are  held 
hxed  at  al^  =  0.14iF^  and  Ki,  =  1.045iFbo-  The  solid  line  in  each  hgure  is  the  predic¬ 
tion  of  Biot  theory  for  the  sintered  glass  beads  with  no  variation  in  the  bulk  modulus  and 
KhQ  =  4.23  X  10®  Pa.  In  (a)  and  (c)  the  sound  speed  ratio  is  dehned  as  the  sound  speed  of 
the  sediment  normalized  by  the  sound  speed  of  the  pore  huid. 
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Eq.  (62)  and  any  larger  values  of  and  would  only  continue  this  trend. 

For  the  sintered  glass  beads,  however,  signihcant  changes  to  the  sound  speed  and  atten¬ 
uation  occur  for  relatively  small  changes  in  the  variance  as  seen  in  Figures  4(a)  and  4(b). 
For  these  values,  the  relation  in  Eq.  (62)  is  easily  satished,  and  the  second-order  perturba¬ 
tion  theory  is  adequate  to  describe  the  mean  field.  In  this  case,  since  the  unperturbed  frame 
modulus  is  so  much  larger  than  for  the  SAX99  sediment,  the  fast  and  slow  compressional 
wavenumbers  are  comparable  in  magnitude  and  the  scattering  from  fast  compressional  waves 
into  fast  compressional  waves  occurs  at  only  slightly  higher  frequencies  than  the  fast  wave  to 
slow  wave  mode  conversion.  While  the  slow  wave  scattering  contribution  to  the  scattering 
levels  out  at  /C2F  ^  10,  the  fast  wave  scattering  attenuation  contribution  causes  the  attenua¬ 
tion  to  continue  increasing  at  a  rate  proportional  to  o;^.  In  the  region  where  the  attenuation 
contributions  from  the  two  scattering  mechanisms  are  comparable,  the  attenuation  increases 
at  a  rate  roughly  proportional  to  uj.  This  is  emphasized  in  Fig.  4(c)  and  4(d)  where  the 
variance  is  held  hxed  and  the  correlation  length  is  varied  from  L  =  0.8  mm  to  L  =  16.6  mm. 
Note  that  due  to  causality,  this  change  in  attenuation  is  also  accompanied  by  a  change  in 
sound  speed.  As  the  slow  wave  scattering  contribution  levels  out,  so,  too,  does  the  sound 
speed.  As  the  contribution  to  attenuation  due  to  the  scattering  of  fast  waves  increases, 
the  sound  speed  begins  to  diverge  as  can  be  seen  in  the  long-short  dashed  in  Fig.  4(c)  for 
L  =  16.6  mm. 

The  second  correlation  function  to  examine  is  the  von  Karman  function, 

C,,  {R)  =  (z/)  {R/iy  K,  {R/L) ,  (65) 


where  is  the  modified  Bessel  function  of  the  third  kind,  T  denotes  the  gamma  function, 
and  u  is  the  Hurst  coefficient  which  is  assumed  to  be  0  <  z/  <  1.  The  three-dimensional 
Fourier  transform  of  this  correlation  function  yields  the  power  spectrum, 

{u  + 1) 


Wtt  (k)  = 


(66) 


7rtr(z/)  (1  +  A;2L2)'^+2 

When  z/  =  1/2,  the  von  Karman  function  reduces  to  the  exponential  correlation  function 
discussed  above. 
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(a) 


(b) 

FIG.  5.  Examples  of  (a)  sound  speed  ratio  and  (b)  attenuation  for  the  sintered  glass  bead 
pack  assuming  a  von  Karman  function  for  the  frame  bulk  modulus  variations.  The  correla¬ 
tion  length  for  each  curve  is  held  hxed  at  L  =  4.1  mm,  the  variance  is  =  0.14iF|,  and 
the  mean  bulk  modulus  is  =  1.045iFfeo-  The  Hurst  coefficient,  z/,  is  z/  =  0.1  (long-dashed 
line),  u  =  0.5  (short-dashed  line),  and  u  =  0.9  (long-short  dashed  line).  In  (a)  the  sound 
speed  ratio  is  dehned  as  the  sound  speed  of  the  sediment  normalized  by  the  sound  speed  of 
the  pore  fluid. 
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Unlike  the  exponential  function,  the  von  Karman  function  does  not  allow  for  an  ana¬ 
lytical  solution  to  the  integral  in  Eq.  (44).  As  a  consequence,  we  will  evaluate  the  integral 
numerically  in  the  subsequent  examples.  Again,  for  the  SAX99  sediment,  the  sound  speed 
and  attenuation  are  only  weakly  affected  for  values  of  and  for  which  the  second-order 
perturbation  theory  applies.  The  variations  in  the  sound  speed  and  attenuation  are  on  the 
same  order  as  those  observed  in  Fig.  3  for  any  value  of  v.  Therefore,  we  will  focus  on  the 
sintered  glass  bead  pack  to  examine  how  varying  v  affects  the  propagating  wave.  Fig.  5 
shows  how  the  sound  speed  and  attenuation  are  affected  as  v  is  varied  from  0.1  to  0.9.  The 
curve  for  v  =  0.5,  the  exponential  function,  is  also  plotted  for  comparison  with  the  previ¬ 
ous  examples.  Increasing  v  above  0.5  has  a  minor  effect  on  the  attenuation,  increasing  the 
scattering  contribution  only  slightly,  while  causing  the  associated  dispersion  to  occur  over 
a  narrower  band  of  frequencies.  The  opposite  occurs  as  v  is  decreased  below  0.5.  In  this 
case,  the  scattering  contribution  to  the  attenuation  decreases  signihcantly  as  does  the  slope 
of  the  attenuation.  The  dispersion  is  also  less  pronounced  as  the  sound  speed  changes  over 
a  much  wider  frequency  band. 


VI.  DISCUSSION 

In  the  introduction,  we  noted  that  for  random  packings  of  granular  materials,  such  as 
sand  grains,  many  simulations  and  experiments  have  observed  the  formation  of  force  chains  in 
the  material.  Along  these  force  chains,  a  subset  of  the  grains  in  the  material  carry  a  majority 
of  the  overburden  pressure  in  the  material.  The  close  relationship  between  the  contact  forces 
between  the  grains  and  the  frame  moduli  suggested  that  the  frame  moduli  should  also  exhibit 
a  random  structure  which  could  lead  to  scattering  losses  for  waves  propagating  through  the 
material. 

However,  for  the  Biot  parameters  measured  for  the  SAX99  sediment,  the  effective  frame 
modulus  is  very  weak  which  is  typical  of  sand  sediments.  In  fact,  a  good  approximation  to 
the  sound  speed  and  attenuation  (neglecting  the  anomalous  high  frequency  behavior  of  the 
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attenuation)  can  be  found  if  the  frame  moduli  are  set  to  zero  as  in  the  Effective  Fluid  Den¬ 
sity  Model  (EDFM)^®.  It  is  not  surprising,  therefore,  that  assuming  a  heterogenous  frame 
modulus  for  this  sediment  should  not  produce  signihcant  contributions  to  the  dispersion  or 
the  attenuation  due  to  mode  conversion  between  the  fast  and  slow  compressional  waves. 
For  both  the  SAX99  and  SAX04  measurements,  as  well  as  the  laboratory  measurements, 
the  attenuation  followed  a  linear  frequency  dependence  and  was  larger  at  high  frequencies 
than  predicted  by  Biot  theory.  The  additional  attenuation  due  to  the  heterogenous  frame 
modulus  is  therefore  insufficient  to  account  for  this  measured  increase  in  attenuation. 

These  scattering  contributions  do,  however,  become  signihcant  for  consolidated  media, 
such  as  the  sintered  glass  bead  example,  where  the  frame  modulus  is  much  larger.  For 
the  sintered  glass  beads,  the  frame  modulus  determined  from  experiments  was  found  to  be 
only  an  order  of  magnitude  smaller  than  modulus  of  grains.  For  the  SAX99  sediment,  the 
frame  bulk  modulus  was  assumed  to  be  at  least  three  orders  of  magnitude  smaller.  For  a 
poroelastic  medium  such  as  the  glass  beads,  where  the  contacts  have  been  fused  together 
to  form  a  fully  consolidated  structure,  the  force  chains  found  in  unconsolidated  media  are 
likely  not  present  since  the  forces  at  the  grain  contacts  have  been  allowed  to  relax  during  the 
melting  process.  The  only  source  of  heterogeneity  in  the  frame  modulus  would  have  to  be 
the  structure  of  the  porous  frame  itself  and  the  local  strength  of  this  kind  of  porous  frame 
is  currently  not  well  understood. 

An  intermediate  case  between  these  two  extremes,  where  both  force  chains  and  a  strong 
frame  modulus  may  be  present,  could  occur  for  a  sediment  under  pressure.  This  could  be  the 
case  for  sand  sediments  tens  of  meters  below  the  sediment/water  interface.  The  overburden 
could  produce  signihcant  contact  forces  between  the  grains  leading  to  a  large  mean  bulk 
modulus  while  the  random  packing  of  the  grains  would  lead  to  heterogeneities  in  the  frame 
modulus. 

To  determine  if  scattering  from  heterogeneities  in  the  sediment  might  lead  to  the  high 
frequency  attenuation  observed  in  the  SAX99  and  in  the  laboratory,  it  remains  to  apply  the 
perturbed  wave  equation  given  in  Eq.  (28)  to  other  heterogeneities  in  the  Biot  parameters 
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listed  in  Table  I.  Of  particular  interest  should  be  those  parameters  for  which  the  predictions 
of  Biot  theory  are  more  sensitive  and  for  which  signihcant  heterogeneities  may  be  expected 
to  occur.  It  is  unlikely  that  the  fluid  parameters  should  vary  signihcantly  and,  while  the 
grains  may  exhibit  some  variation  in  their  modulus  and  density,  sorting  due  to  wave  action 
may  lead  to  a  suppression  in  the  variation  between  the  grains  themselves.  This  suggests 
that  those  parameters  which  relate  directly  to  the  pore  structure  itself,  specihcally  the 
permeability,  tortuosity,  and  the  porosity,  may  exhibit  spatial  variations  which  could  lead  to 
signihcant  scattering.  Work  is  currently  underway  to  examine  the  predictions  of  the  general 
perturbation  theory  given  in  Eq.  (28)  when  the  porosity  of  the  sediment  is  heterogenous^^. 


VII.  CONCLUSIONS 

In  an  effort  to  understand  the  role  of  heterogeneities  in  the  propagation  of  high-frequency 
sound  through  an  ocean  sediment,  we  have  applied  perturbation  theory  to  Biot’s  poroelastic 
equations.  The  resulting  wave  equation,  while  capable  of  modeling  heterogeneities  in  any  of 
the  Biot  parameters,  was  applied  to  the  case  of  a  randomly  varying  frame  bulk  modulus.  The 
presence  of  the  heterogeneities  in  the  bulk  modulus  affects  both  fast  and  slow  compressional 
wave  propagation,  but  does  not  affect  the  propagation  of  shear  waves.  At  low  frequencies,  the 
resulting  mean  wavenumber  for  the  fast  compressional  wave  reduces  to  the  result  presented 
previously  by  Miiller  and  Gurevich  for  seismic  propagation^®. 

The  effects  of  the  frame  modulus  heterogeneities  can  be  incorporated  into  a  set  of  effec¬ 
tive  Biot  coefficients  each  of  which  consists  of  the  mean  coefficient,  a  frequency-independent 
term  that  is  a  function  of  the  frame  modulus  variance,  and  a  frequency-dependent  term 
that  is  a  function  of  the  frame  modulus  covariance  within  the  medium.  The  frequency- 
independent  term  accounts  for  the  softening  of  the  medium  due  to  the  variations  in  the 
frame  modulus,  while  the  frequency-dependent  accounts  for  scattering  from  the  hetero¬ 
geneities.  For  a  propagating  fast  wave,  when  k2L  ^  1,  mode  conversion  between  the  fast 
and  slow  waves  scatters  energy  from  the  mean  held  increasing  the  attenuation  of  the  mean 


field.  At  higher  frequencies,  energy  is  also  lost  due  to  scattering  from  the  coherent  fast  wave 
into  incoherent  fast  waves.  While  the  attenuation  due  to  mode  conversion  levels  off  near 
k2L  ^  10,  the  attenuation  due  to  fast  wave  scattering  goes  as 

The  perturbation  theory  accounts  for  scattering  up  to  second  order  in  the  modulus 
variations  and  was  shown  to  be  applicable  to  small  variations  in  the  bulk  modulus.  For 
weakly  consolidated  granular  media,  such  as  sand  sediments,  the  mean  frame  modulus  is 
typically  much  smaller  than  the  grain  or  fluid  modulus  and  the  role  of  the  variations  in  the 
frame  modulus  were  shown  to  only  weakly  affect  the  sound  propagation.  For  consolidated 
media,  such  as  the  sintered  glass  bead  pack  discussed  here,  weak  variations  in  the  frame 
modulus  were  shown  to  have  a  significant  effect  on  sound  propagation.  Future  work  will 
explore  the  role  of  heterogeneities  in  parameters  such  as  the  porosity,  which  may  have  a 
more  signihcant  affect  on  propagation  through  weakly  consolidated  poroelastic  media. 


APPENDIX  A:  BIOT  GREEN’S  TENSOR 

In  order  to  apply  perturbation  theory  to  solve  (16)  with  varying  frame  bulk  and  shear 
moduli,  we  need  to  know  the  Green’s  function  which  can  be  used  to  solve  the  inhomogenous 
poroelastic  equation, 

[Ko]V(V-U)-[po]Vx(VxU)  + 

[po]o;2U=-S(r).  (Al) 

The  inhomogenous  equation  with  a  point  source, 

[Ko]  V  (V  ■  0  (R))  -  [po]  V  X  (V  X  0  (R))  + 

[Po]o;20(R)  =  -<5(R)[|]  a,  (A2) 

where  [I]  is  the  2x2  identity  matrix,  3  is  the  identity  tensor,  and  R  =  |r  —  r'|,  has  the 


29 


Green’s  tensor  solution 


©(/?)  =5^[G„]  VVgniR) 

n=l 

-  [G,]  V  X  V'  X  (g,  (R)  3) 

-[po]-^u;-^5{R)3,  (A3) 


where  n  =  1  and  2  is  associated  with  the  fast  and  slow  wave,  respectively. 

s  is  associated 

with  the  shear  wave. 

,,  ,  [K,]-^-[Po]-^u;-^kl 

-  P  -  P 

/il  /i2 

(A4) 

[Ko]-^-[po]-^u;-^kl 
^2  ^1 

(AS) 

lG.]  =  h^+|Po]-'i*'X 

(A6) 

^ihmR 

gm  {R)  =  where  m  =  1,  2,  or  s. 

(A7) 
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-P*  Pf 
-Pf  P 


ip)  -  pp’) ' 


(A9) 


and 


0  0 


(AlO) 


L°  aJ 

Note  that  the  notation  for  the  operators  in  (A3)  differs  from  the  notation  of  Morse  and 
Feshbach  who  write  them  as  Vg'n  {R)  V'  and  V  x  gn  {R)  x  Vh  Morse  and  Feshbach  illustrate 
the  meaning  of  W'gn  {R)  by  operating  with  W'gn  {R)  on  some  vector  f  (r')  which  yields, 


fa:  (r')  V  {dg/dx’)  +  fy  (r')V  {dg/dy') 

+Mr')V{dg/dz'). 


(All) 
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The  other  operator  can  be  understood  by  writing  it  as 


V  X  V'  X  (ff,  (R)  3)  =  VV'g,  (R)  +  3V\  (R) . 


(A12) 


Also  note  that 


1  dgn  _ 


V7V7^  (  TD\  I  ^00  ^  r-i 


(A13) 


where  R  is  the  unit  vector  in  the  direction  of  R. 

If  the  source  term  in  (Al)  can  be  separated  into  longitudinal  and  tangential  components, 


S  (r)  =  V<h5  +  V  x^s, 


(A14) 


the  resulting  held  can  be  determined  using  the  Green’s  tensor  as 
2 


U  (r)  =  ^  kl  [G„]  V  f  Qn  (R)  (r')  dh'  + 

n=l  d 

kl[Gs\Vx  j  g,  (R)  (rOdV  + 
V  X  J  S{R)^s  (r')  dhk 


Two  relations  that  are  used  to  simplify  the  perturbation  solution  are 


Z  [GJ  =  [Ko]-' 


n=l 


and 


EiGJ  =  [p»i'‘ 


-1  -2 

OJ  . 


(A15) 


(A16) 


(A17) 


n=l 


APPENDIX  B:  MEDIUM  SOFTENING  DUE  TO  MODULUS  VARIATIONS 


To  understand  the  origins  of  the  frequency  independent  term  in  the  effective  Biot  co¬ 
efficient  matrix,  Eq.  (38),  we  consider  the  mean  response  of  a  poroelastic  medium  under 
static  compression.  For  the  heterogenous  poroelastic  medium  we  are  considering,  we  will 
take  the  harmonic  mean  of  the  Biot  coefficient  matrix  to  determine  the  effective  response  of 
the  medium, 

|K.„]-'  =  (|K(r)|-'>.  (Bl) 
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To  see  that  this  is  equivalent  to  the  self-consistent  method  for  estimating  the  elastic  constants 
in  inhomogenous  material  discussed  by  Berryman^^,  consider  what  happens  if  we  replace  the 
fluid  in  the  poroelastic  medium  with  a  vacuum.  In  this  case,  the  Biot  coefficients  reduce 
to  C  =  M  =  0  and  H  (r)  =  Kf,  (r)  -|-  4/ife/3,  pf  =  0,  and  the  Biot  equations  reduce  to  the 
equations  for  an  elastic  medium.  The  harmonic  mean  of  H  (r)  becomes 


H. 


eff 


H(r] 


1 


KJr)  +  ^r 


(B2) 


‘■bi-y  '  3 

which  is  equivalent  to  Eq.  (45)  in  reference 

With  variations  in  the  frame  bulk  modulus  in  the  form  of  K/,  (r)  =  Ki,q  +  SKi,  (r),  we 
can  approximate  [Keff]  by  expanding  [K  (r)]~^  in  a  power  series  about  Ki,q, 


[K(r)]-^  =  [Ko]-^  + 


-1  , 


-1 


dK, 


-5Kh  (r) 


(B3) 


where  we  have  retained  only  those  terms  up  to  second  order  in 

The  partial  derivatives  of  the  inverse  Biot  coefficient  matrix  in  Eq.  (B3)  can  be  related 
to  the  partial  derivatives  of  the  Biot  Coefficient  matrix  by 


OKt, 


=  -[Ko 


OKt, 


[K( 


T-l 


(B4) 


[N  _  ri,  1-1  <9^  [Ko]^  1-1 


OKI 


=  -[Kc 


dK. 


[Ko 


(B5) 


aiit  ‘  SAl 

Substituting  Eq.  (B4)  and  (B5)  into  Eq.  (B3)  and  taking  the  mean  of  the  expansion  gives 


1-1  ri/  1-1  9  ri/  1-1  ^  riy  1-1  _2 


[Keff]  -  [Ko]  -  -  [Ko] 


2  OKi 


[Ko]  (Tl 


bb 


ru'  1—1  ^  [Kq]  ]-!  *9  [Kq]  .|_i  ^2 
[Ko]  [Ko]  [Ko]  CTj 


bb- 


(B6) 


dKk  ^  dK, 

To  hnally  arrive  at  an  expression  for  [Keff],  we  need  to  determine  the  inverse  of  Eq. 
(B6).  To  do  this  we  will  expand  the  inverse  of  the  right  side  of  Eq.  (B6)  as  a  function  of  the 
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variance  about  =  0.  Upon  doing  this  and  using  the  relationship  between  the  derivative 
of  a  matrix  and  the  derivative  of  its  inverse  illustrated  in  Eq.  (B4),  we  arrive  at  the  hnal 
approximation  for  the  effective  Biot  coefficients, 
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This  expression  is  identical  to  the  frequency  independent  term  in  Eq.  (38)  determined  from 
perturbation  theory. 
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TABLE  I.  Parameters  used  in  calculating  Biot  results  unless  otherwise  stated.  The  values  in 
the  column  of  parameters  labeled  “SAX99  Sand”  correspond  to  the  values  in  Table  I  of  Ref. 
6  for  the  best  £t  to  the  sound  speed  and  attenuation  measured  during  SAX99.  Note  that 
for  the  frame  moduli,  the  imaginary  parts  have  been  removed  and  the  moduli  are  assumed 
to  be  purely  real.  The  values  in  the  column  of  parameters  labeled  “Sintered  Glass  Beads” 
were  determined  from  the  parameters  listed  for  the  sintered  glass  bead  sample  B  in  Table 
III  of  Ref.  30. 


Parameter 

SAX99  Sand 

Sintered  Glass  Beads 

Porosity  {{3) 

0.385 

0.37 

Mass  density  of  sand  grains  {pg) 

2690  kg/m^ 

2530  kg/m^ 

Mass  density  of  water  (pj) 

1023  kg/m^ 

1023  kg/m^ 

Bulk  modulus  of  sand  grains  {Kg) 

3.2  X  10^°  Pa 

4.3  X  10^°  Pa 

Bulk  modulus  of  frame  (iL;,) 

4.36  X  10^  Pa 

4.23  X  10^  Pa 

Shear  modulus  of  frame  (pb) 

2.92  X  10^  Pa 

7.94  X  10^  Pa 

Bulk  modulus  of  water  {Kj) 

2.395  X  10*^  Pa 

2.395  X  10^  Pa 

Viscosity  (p) 

0.00105  kg/m-s 

0.00105  kg/m-s 

Permeability  {k) 

2.5  X  10-1^  m2 

4.43  X  10-1^  m2 

Tortuosity  (a) 

1.35 

1.60 
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FIG.  1  (a)  Fast  compressional  sound  speed  and  (b)  attenuation  determined  from 

Eq.  (41)  with  [Keff]  given  by  Eq.  (54)  as  a  function  of  the  variance, 
of  the  moduli  variations  normalized  by  K"^.  Both  the  sound  speed  and  the 
attenuation  have  been  normalized  by  their  respective  values  when  abb  =  0. 

These  results  were  calculated  for  /  =  1  kHz  using  the  SAX99  sediment 
parameters  in  Table  I  while  the  value  of  the  mean  bulk  modulus  was  increased 
from  4.36  x  10^  Pa  for  the  solid  line,  to  4.36  x  10®  Pa  for  the  short  dashed 
line,  and  to  2.18  x  10®  Pa  for  the  long  dashed  line .  14 


FIG.  2  Values  of  al^/ for  which  [Keff]  remains  equal  to  [Kq]  for  the  SAX99  sediment 
parameters  (solid  line)  and  the  sintered  glass  bead  parameters  (dash-dot  line) 
as  Kb  is  increased  by  three  orders  of  magnitude .  16 

FIG.  3  Examples  of  (a)  sound  speed  ratio  and  (b)  attenuation  for  the  SAX99  sedi¬ 
ment  assuming  an  exponential  correlation  for  the  frame  bulk  modulus  varia¬ 
tions.  The  correlation  length  is  held  hxed  at  L  =  1  mm  while  the  variance  and 
the  mean  bulk  modulus  are  increased.  The  variance  is  shown  in  the  figures 
and  the  mean  bulk  modulus  is  Kb  =  l.lSKbo  (long-dashed  line),  Kb  =  lAlKbo 
(short-dashed  line),  and  Kb  =  2.041^60  (long-short  dashed  line).  In  (a)  the 
sound  speed  ratio  is  dehned  as  the  sound  speed  of  the  sediment  normalized 
by  the  sound  speed  of  the  pore  fluid.  In  (b)  the  lower  set  of  curves  is  the 
contribution  of  the  scattering  from  the  fast  compressional  waves  into  the  slow 
compressional  waves  to  the  attenuation  of  the  mean  fast  compressional  wave. 

The  solid  line  in  each  hgure  is  the  prediction  of  Biot  theory  with  no  variation 
in  the  bulk  modulus  and  Kb  =  Kbo .  21 
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FIG.  4  Examples  of  sound  speed  ratio  and  attenuation  for  the  sintered  glass  beads 
assuming  an  exponential  correlation  for  the  frame  bulk  modulus  variations. 

In  (a)  and  (b),  the  correlation  length  is  held  fixed  at  L  =  1.6  mm  while  the 
variance  and  the  mean  bulk  modulus  are  increased.  The  variance  is  shown  in 
the  hgures  and  the  mean  bulk  modulus  is  =  l.OOQi^w  (long-dashed  line), 

Kij  =  1.021iFfeo  (short-dashed  line),  and  =  1.045iFbo  (long-short  dashed 
line).  In  (c)  and  (d),  the  variance  and  mean  bulk  modulus  are  held  hxed 
at  =  0.14iF^  and  =  l.OdSiF^o-  The  solid  line  in  each  hgure  is  the 
prediction  of  Biot  theory  for  the  sintered  glass  beads  with  no  variation  in  the 
bulk  modulus  and  =  4.23  x  10®  Pa.  In  (a)  and  (c)  the  sound  speed  ratio 
is  defined  as  the  sound  speed  of  the  sediment  normalized  by  the  sound  speed 

of  the  pore  fluid .  23 

FIG.  5  Examples  of  (a)  sound  speed  ratio  and  (b)  attenuation  for  the  sintered  glass 
bead  pack  assuming  a  von  Karman  function  for  the  frame  bulk  modulus 
variations.  The  correlation  length  for  each  curve  is  held  fixed  at  L  =  4.1  mm, 
the  variance  is  =  0.14iF|,  and  the  mean  bulk  modulus  is  =  1.045i^w 
The  Hurst  coefficient,  z/,  is  z/  =  0.1  (long-dashed  line),  u  =  0.5  (short-dashed 
line),  and  u  =  0.9  (long-short  dashed  line).  In  (a)  the  sound  speed  ratio  is 
dehned  as  the  sound  speed  of  the  sediment  normalized  by  the  sound  speed 
of  the  pore  fluid .  25 
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